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A simple model of a degenerate two-dimensional Bose liquid interacting with a fluctuating gauge 
field is investigated as a possible candidate to describe the charge degree of freedom in the normal 
state of the cuprate superconductors. We show that the fluctuating gauge field efficiently destroys 
superfluidity even in the Bose degenerate regime. We discuss the nature of the resulting normal 
state in terms of the geometric properties of the imaginary-time paths of the bosons. We will also 
present numerical results on the transport properties and the density correlations in the system. We 
find a transport scattering rate of %/tu ~ 2ksT, consistent with the experiments on the cuprates 
in the normal state. We also find that the density correlations of our model resemble the charge 
correlations of the t-J model. 
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We study the low-temperature behavior of repulsive 
bosons in a spatially fluctuating gauge field in two di- 
mensions. This is motivated by the gauge theories of the 
t-J model for the cuprate superconductors, where low- 
energy charge excitations are described by bosonic de- 
grees of freedom. The internal gauge field of this model 
suppresses superfluidity in the Bose liquid, even below 
the Bose degeneracy temperature when there is signif- 
icant exchange among the bosons. We can study the 
imaginary-time trajectories of the bosons in the path- 
integral representation of this model. We see that the 
boson world-lines retrace themselves in the presence of 
strong gauge fluctuations, giving rise to interesting dy- 
namics in this degenerate but metallic Bose liquid. 

We have studied this metallic state using quantum 
Monte Carlo techniques. We find that this model does 
indeed capture some of the long-wavelength charge prop- 
erties which are common to the cuprate superconduc- 
tors. This includes a linear temperature dependence of 
the transport scattering rate l/r tr , as deduced from a 
Drude-like optical conductivity from our model. This is 
consistent with experimental data on the cuprate super- 
conductors near optimal doping. We also find that the 
density excitations in our model are qualitatively similar 
to those in the full t-J model, by comparing our results 
with diagonalization results in the literature. A brief ac- 
count of this work has already appeared!!. 



I. MOTIVATION 

The normal metallic state of the superconducting 
cuprates displays many non-Fermi-liquid properties. For 
instance, the in-plane resistivity of La2- a: Sr a ;Cu04 has a 
power-law temperature dependence of the form p cx T a 
where *at increases from 1 to 1.5 with increasing hole 
dopingp. In particular, near optimal doping, the resistiv- 
ity is linear in temperature up to 1000K. This linear-T 
dependence is found in many of the cuprate superconduc- 
tors with similar values of dp/dT (1.2/xQcm/K ± 2O%)0. 



This should be contrasted with the quadratic tempera- 
ture dependence of Fermi-liquid theory. Similarly, the 
transport relaxation rate appears to be universal among 
optimally-doped compounds: l/r tr ~ 2k^T (from a two- 
component-model analysis of the .optical conductivity in 
YBCOq, LSCOB, Bi2212B, Bi2201Q). Transport in a mag- 
netic field is also anomalous. The Hall coefficient in- 
dicates the existence of hole-like carriers in the doping 
range where superconductivity occurs. The Hall coeffi- 
cient i?u increases with decreasing temperature, but it re- 
mains smaller than the classical value of 1 /rihec for a hole 
density of rit for a wide range of temperatures down to 
the superconducting transition. These compounds also 
have a small positive magnetoresistance with a tempera- 
ture dependenceu different from conventional theory us- 
ing Ttr- 

The transport properties of these compounds appear 
to have common features in spite of considerable differ- 
ences in the transition temperature and spin fluctuation 
properties among these compounds. This indicates that 
a common mechanism is responsible for the scattering of 
charge carriers in these materials. One might hope that 
this scattering mechanism can be understood in terms of 
a low-energy theory with a minimum number of micro- 
scopic parameters. In this paper, we study a Bose liquid 
in a fluctuating gauge field as a possible candidate for 
such an effective theory. 

The anomalous transport behavior, together with 
other unusual features such as temperature-dependent 
magnetic susceptibility and non-Korringa behavior of the 
nuclear magnetic relaxation time, leads to the conclu- 
sion that the metallic state of the cuprates cannot be 
described in a simple Fermi-liquid scenario. It has been 
postulated that "spin-charge separation" is responsible 
for these anomalies!. For instance, such a scenario might 
reconcile the apparent low density and hole-like charac- 
ter of the charge carriers with the observation of a large, 
electron-like Fermi surface in photoemission. Numerical 
studies of the t-J model, which is believed to be a low- 
energy model of the cuprates, also provide some support 
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for spin-charge separationBiLil, such as different energy 
scales for the spin and charge excitations, and the sup- 
pression of 2/ci?-scattering in the charge spectrum. 

A model of spin-charge separation is a gauge theory 
where neutral spin-half fcrmions ( "spinons" ) and charge- 
e bosops ("holons") interact via an internal U(l) gauge 
fieldlI3't± Physically, the transverse part of_thc gauge 
field is related to "spin chirality" fluctuationsEl In this 
picture, the charge properties of the system should be 
dominated by the behavior of the holons. We will study 
the holon subsystem in this paper, treating the spinon 
subsystem simply as a medium through which the gauge 
field propagates. To be more precise, we study a model 
of bosons with on-site repulsion in the presence of a spa- 
tially fluctuating magnetic field with short-range corre- 
lations. The repulsive interaction is necessary for the 
stability of the system, which means that one cannot 
treat this problem perturbatively starting from an ideal 
Bose gas. Previous studiesEJlia have implicitly studied 
the non-degenerate regime of low density or high tem- 
perature, whereas the regime relevant to the cuprates 
is the degenerate regime where the thermal deBroglie 
wavelength of the bosons is greater than the mean parti- 
cle spacing. A concern from earlier studies of the gauge 
model is that degenerate bosons would have strong dia- 
magnetic response to the internal gauge field and hence 
effectively Bose-condense at a relatively high tempera- 
ture {k B T B E ~ imht ~ I000K). This would in fact re- 
store Fermi-liquid behavior to the system. We shall show 
here that gauge fluctuations suppress this diamagnetic 
response and the bosons remain normal without strong 
diamagnetism at all finite temperatures. Furthermore, 
our numerical results indicate that the resistivity of this 
Bose metallic phase has a linear temperature dependence 
which is consistent with experiments. 

Besides the possible relevance to the transport in the 
cuprate superconductors, the model we consider is of in- 
trinsic theoretical interest. The model is a Bose version of 
the problem of a quantum particle in a random magnetic 
flux, which has received considerable attention in recent 
years. It is also related to frustrated spin systems and 
vortex glasses. However, since we deal exclusively with 
annealed averaging in this paper (see later), we cannot 
draw any direct conclusions about these problems with 
quenched disorder. 

The rest of the paper is organized as follows. In section 
II, we review the connection between the gauge theory 
of the t-J model and our boson model. In section III, 
we discuss the path-integral formalism which provides a 
convenient framework to visualize physical processes in 
terms of the imaginary-time paths of the bosons. In sec- 
tion IV, we look at the effects of the gauge field on the 
world-line geometry of the bosons. We will see that the 
partition function of the system is dominated by self- 
retracing world-line configurations. We will also argue 
that superfluidity is destroyed by the fluctuating gauge 
field, giving rise to a degenerate Bose metal. In the sub- 
sequent sections, we present the results of a quantum 



Monte Carlo study of this metallic phase. We will dis- 
cuss the transport properties and the density correlations 
in this boson model. 



II. A BOSON GAUGE MODEL 

In this section, we provide the motivation for studying 
an effective boson model from the gauge theory of the 
t-J model, which describes the motion of vacancies in a 
doped Mott insulator: 

n = -t ]T (cU- CT + h.c.) + • s j (!) 

with the constraint of no double occupancy. Experimen- 
tally, J ~ 1500K and t Q /J ~ 3. 

The constraint of no double occupancy allows us to 
write the creation of a physical hole in terms of the cre- 
ation of a charged hard-core boson (holon) and the anni- 
hilation of a spin-half fermion (spinon): Ci a — fi a b\. In 
terms of these slave bosons and fermions, the Hamilto- 
nian of the t-J model can be written as: 

n = -*° E (flMfj- + h - c -) • -'E 8 -- s < 

+ J2 a oi(flfi- + bth-l) (2) 

i 

where S, = f} a cr a pfi0> The aoi-field is a Lagrange multi- 
plier enforcing the local occupancy constraint, and acts as 
a fluctuating scalar potential for the spinons and holons. 

Among the mean field theories proposed to decouple 
the quartic terms in Eq. (Q), a candidate for the normal 
state near optimal doping is the the uniform resonating- 

valence-bond (RVB) ansatz: E CT (/t/jv) = ^ 3 ■ This 
incorporates short-range antifcrromagnetic correlations 
without any long-range Neel order. The Lagrangian of 
this RVB phase can be written as: 

C = ^2 fia(dr - MF + aot)/i<7 + X! b *( dr ~ ^ B + a Qi) b i 

1,(7 i 

m 

-t £5> to «&;6 i + h.c.) (3) 

(ij) 

The vector potential arises from the fluctuations 
in the phase of the RVB order parameter. Longitudi- 
nal fluctuations of the gauge field ay do not affect the 
Lagrangian due to an internal U(l) gauge symmetry: 
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We will therefore work in a fixed gauge, such as the 
Coulomb gauge, and consider only the fluctuations in the 
transverse part of the gauge field a.^ . In other words, we 
will consider only fluctuations in the internal magnetic 
and electric fields which are gauge- invariant quantities. 

Since we are interested in the charge degrees of free- 
dom, we wish to consider an effective theory with bosons 
only, and regard the spinon fluid as a medium through 
which the gauge field propagates. The gauge field has 
no dynamics in vacuo. The response of the spinon fluid 
to the gauge field is responsible for the dynamics of the 
gauge field as seen by the holons. More specifically, we 
can obtain the Gaussian fluctuations of the a-fields by 
treating the spinon response in the random-phase ap- 
proximation. The effective gauge-field propagator is: 

Sg = n ° (k,^rO a o(k,^n)ao(k,cj n ) 



|^ ^] n 1 (k,w„)a^(k,w n )a i (k,w n ), (5) 



2/3L 



where /3 = 1/T, u) n = 2 / nnT 1 L is the linear size of the 
system, and a± is the transverse part of the gauge field. 
(We use units where distance is measured in terms of the 
lattice spacing and fee = h = e = 1.) Here, for small 
k and u n , II 00 ~ pp, the spinon density of states at the 
Fermi level. This describes the Thomas-Fermi screen- 
ing of internal electric fields by the fermions. The ef- 
fective interaction mediated by the screened ao-field is a 

— 1/2 

repulsion between the bosons (of range oc p F ), con- 
sistent with the original hard-core requirement for the 
bosons. We will model this with an on-site repulsion en- 
ergy, U. On the other hand, the magnetic fields due to 
fluctuations in are not effectively screened out by the 
fermionsll3. The gauge-field fluctuations as experienced 
by the holons are therefore strong. More specifically, 
the Gaussian fluctuations have the correlation function 
-D(k, u> n ) = (dj_(k, cj n )a_i_(k, w n )), given (in the contin- 
uum limit) by: 



1 



1 



II_L(k,w n ) ~/\u; n \/k + xk 2 



(6) 



where x is t ne orbital susceptibility of the spinon fluid 
and 7 is a Landau damping coefficient. These gauge- 
field fluctuations cause profuse forward scattering of the 
bosons. We believe that this is the dominant scattering 
mechanism in this problem. Since it is overdamped at 
long wavelengths with a relaxation rate which diverges 
as 1/fc 3 , we will ignore the slow relaxation and work in a 
"quasistatic" limit for the gauge fields: 



D(k, u n ) -»■ D(k, uj n = 0) 5 n>0 = 5 nfi /xk 2 



(7) 



(On a square lattice, k 2 is replaced by I — (cosk x — 
cosfc y )/2.) This quasistatic approximation is justified 
when the gauge field relaxes on a time scale longer than 
1/T. Since the typical scattering wavevector of interest is 



the inverse deBroglie wavelength of the bosons, the rele- 
vant relaxation time scales as 1/k 3 ~ 1/T 3 / 2 . One might 
therefore expected that this approximation is valid at a 
sufficiently low temperature. 

One might object that arguments above are based on 
a weak-coupling theory of the response of the spinons to 
the gauge fields. However, we believe that the essential 
features remain correct in general, namely a separation 
of times scales between the relaxation of the gauge fields 
and the boson dynamics, as well as the magnitude of the 
gauge fluctuations being controlled by the spinon dia- 
magnetic susceptibility x- 

The gauge-field correlator (Q) corresponds to a spa- 
tially uncorrelated flux distribution with the correlation 
function: 



T 

($ r $ r ,) = -5 r y 

X 



(8) 



where $ r = ($ /2vr) ^ n a,ij (oriented sum around the 
links of plaquette r) is the flux through plaquette r. 
($0 = hc/e is the flux quantum.) Since we are treat- 
ing the thermodynamics for the gauge field classically, 
we have a thermal factor of T in Eq. (|J) for the flux 
variance ( < I )2 ). Given that the fermion orbital suscepti- 
bility is roughly constant at low temperatures, we might 
expect the flux variance to have a linear temperature 
dependence. However, a lattice calculation by Hlubina 
et alEJl has indicated that the Gaussian fluctuations are 
sufficiently strong that the flux through a plaquette is 
of the order of the flux quantum <& : ($ 2 ) : / 2 > 0.5<!>o 
down to a temperature of 0.4J. Since the experimen- 
tal superconducting T c is of the order of 0.1J, we ex- 
pect that this regime of strong random flux is relevant 
to the normal state of the cuprates until one approaches 
the superconducting transition. In this regime, the pre- 
cise value of ($ 2 ) should not affect the behavior of the 
bosons, and we will focus on a large and temperature- 
independent flux variance when we study the transport 
and correlation functions of our boson system. 

Another factor leading to the reduction of the flux 
variance at low temperature is one that has not been 
discussed so far, namely that the magnitude of the 
gauge field should also be affected by the diamagnetic 
response of the holons as well as the spinons, i.e., ($ 2 ) = 
T/(x S pinon(T) + Xhoion(T)). The holon contribution dom- 
inates near an instability to Bose condensation where 
Xhoion diverges and the bosons develop a Meissner re- 
sponse to expel the gauge field from the system alto- 
gether. However, we will see in this paper that Bose 
condensation and the holon diamagnetism are strongly 
suppressed even below the boson degeneracy tempera- 
ture. Therefore, in a wide range of temperatures above 
the superconducting T c , we are justified in neglecting this 
feedback effect of the holons on the magnitude of the 
gauge field fluctuations. 

We can now define more precisely the effective model 
which we study in the rest of the paper. It is a model of 
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lattice bosons interacting with a quasistatic gauge field, 
described by effective action 5 = 5b + 5q ■ 



Sb = b*(d T - MB )^ - H b (t)^J dr 



k 

2{¥) 



(9) 
(10) 



with the boson Hamiltonian 



H B = -i^y^bfo +h.c.) + - 1), (11 



where t = tot; ~ to, L is the linear size in units of lattice 
spacing and U ^S> t. Note that, on performing the average 
over the gauge field, we average over static configurations 
only, i.e., a(k,w„ ^ 0) = 0. 

We cannot say that we have rigorously derived above 
effective action from the slave-boson mean field theory 
of the t-J model. Many approximations have been intro- 
duced to obtain this simple model with few adjustable 
parameters. For example, we have neglected the tem- 
perature dependence of the RVB order parameter-/ and 
also the gauge-field correlations of higher ordeiS. We 
take the point of view that we are studying a "minimal" 
low-energy theory which hopefully captures many of the 
generic features of more complicated models. 



III. PATH INTEGRAL REPRESENTATION 

It is convenient to study our boson model in a first- 
quantized formulation. The partition function Z for a 
system with N bosons in the canonical ensemble can be 
written in terms of a Feynman path integralcJ over the 
boson trajectories {x a (r)} (a = 1, . . . , N): 



E 



c(0)=P(x(/3)) 

D[xi, ...,Xjv] x 



(12) 



where S B is the action for bosons in the absence of mag- 
netic fields: 



ch 



b*d T bi 



(13) 



where iJg is given by (11) with =0. In this sec- 



tion, we will discuss the model in the continuum limit 
for notational convenience. In the continuum, one would 
have: 



dr 



' N 

Em . n 

a=l 



+ VC/5(x Q (r)-x 7 (r)) 



a>7 



(14) 



Particle identity is taken into account by perform- 
ing the path integral over all trajectories where the 
set of final boson coordinates at {xi (/?),... ,xjv(/3)} 
is some permutation of the initial boson coordinates 
{xi(0), . . . , xjy(O)}. Any such permutations can be bro- 
ken down to cycles. Each cycle forms a closed loop when 
the imaginary-time trajectories (world lines) of a many- 
boson configuration are projected onto real space. At 
high temperatures, cycles of length 1 dominate the par- 
tition function and the system is in a non-degenerate 
classical regime. At temperatures below the degener- 
acy temperature of the bosons, particles can travel large 
distances in the imaginary time, forming many ring ex- 
changes (see Fig. [l]). 




FIG. 1. A schematic configuration for 6 bosons after pro- 
jecting the imaginary-time paths onto the xy-p\&ne. There 
are a total of 3 cycles: 1 cycle of one particle, 1 cycle of two 
particles, and 1 cycle of three particles. Solid circles denote 
particle positions at r = and /3. 

In this formulation, we may integrate out the Gaus- 
sian fluctuations of the gauge field in ([l2]). Thus, we 
arrive at a boson-only effective theory which we study 
numerically in this work. The system is described by 
the partition function Z = J 2?xe~ Sotf where the effective 
action is given by: 



S e s = 5* B + 52 



(15) 



with 



S 2 =iV/ / D{* a {T)-* a ,{T'))± a -± a ,dTdT'. (16) 
act' J ° J ° 

where D(x) = (1//3L 2 ) £ k #o D ( k > 0)e" lk x . Note that 
the k = contribution has been excluded in the sum 
over k, corresponding to a gauge choice where the k = 
part of a is zero. This is one way to fix the remaining 
degree of gauge freedom which is not determined by the 
condition of V • a = 0. If we consider a system with 
periodic boundary conditions in space, another scheme 
would be to fix the line integral of the gauge field around 
a specified path which wraps around the boundary. How- 
ever, the latter scheme is inconvenient for our purposes 
because it breaks translational invariance explicitly. 

The current interaction .D(x) mediated by the gauge 
field is logarithmic at large distances, and is attractive 
between opposite currents. Due to the quasistatic na- 
ture of the gauge fields, the interaction is also infinitely 
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retarded in time. We will see in the next section that 
this encourages world lines to retrace themselves, with 
important consequences for the boson dynamics. 

Before proceeding to discuss the physical consequences 
of the current interaction 62 , some remarks about our av- 
eraging procedure for the gauge fields are in order. We 
have performed an "annealed" average over the gauge 
fields, rather than a "quenched" average. Annealed av- 
eraging is necessary in our case because our gauge field 
a is an internal thermodynamic variable. Formally, we 
evaluate observables (O) as: 



1 



VxOe 



-Sett _ 



/PxPaP[a]0e' 5g ''/ adx 



(17) 



where P[a] = Af~ 1 5(V • a)e~ SG [ a l is the probability dis- 
tribution for the gauge field, and J\f is a suitable normal- 
ization factor. This is different from quenched averaging 
which would be appropriate if we dealt with a system 
with frozen impurities, such as a vortex glass. Quenched 
averaging requires the evaluation of: 



VaP[a] 



JVxO 



,-Ss-t J adx 



/Pxe" 5 B- l / a dx . 



(18) 



The differences between quenched and annealed averag- 
ing from the point of view of perturbatio n .. (d iagram- 
matic) theory has been addressed elsewheretlrEJ. 

From the point of view of the path integral Monte 
Carlo method, our ability to perform the annealed av- 
eraging means that we would not have to perform ex- 
tensive averages over different frozen realizations of the 
random flux. Moreover, note that the effective action 
( |l5| ) is manifestly real, and so we avoid the sign problem 
which occurs numerically when performing a quenched 
average over the gauge fields. We have studied boson 
densities between rib = 1/4 and 1/6. We choose an on- 
site interaction strength U > At. We follow the Moiate 
Carlo methods of Ceperley and Pollockt3 and Trivedica 
Each Monte Carlo step involves the reconstruction of the 
world lines, {x Q (r)}, for all N particles using the ideal 
boson propagator in a short interval in imaginary time. 
The on-site interaction and the current interaction S% 
are taken into account using Metropolis tests. To en- 
sure quantum exchange, we may insist that each accepted 
configuration differs from the previous one by a pair ex- 
change. This can be incorporated, without loss of de- 
tailed balance, as a Metropolis test. We refer readers to 
the original referencesEJcj for further details. (In evalu- 
ating the gauge field contribution S2 , we have also made 
use of a geometrical interpretation of S2 which we discuss 
in the next section.) In the discretization of the imag- 
inary time, we have used a small At = j3/M < 0.1/t, 
so as to minimize the systematic error and to allow the 
reliable use of maximum entropy techniques to perform 
analytic continuation on our imaginary-time data to ob- 
tain the dynamical quantities of interest. This sets the 



lowest accessible temperature to T ~ O.lf for lattice sizes 
considered here. For studies on dynamic response to be 
discussed later, we have restricted ourselves to lattices of 
sizes up to 6 x 6, due to the need to obtain imaginary 
time correlation functions to a high accuracy. For the 
calculation of static properties, we have studied lattices 
as large as 10 x 10. 

To summarize, we have obtained an effective theory of 
bosons with current interactions which are long-ranged 
in space and time. This model can be studied using path 
integral Monte Carlo methods. In the next section, we 
will discuss how these interactions affect the geometry of 
the boson world lines and hence the physical properties 
of the system. 



IV. EFFECT OF GAUGE FIELDS ON WORLD 
LINE GEOMETRY 

A. "Brinkman-Rice bosons" 

In this section, we will discuss how the current inter- 
action 5*2 mediated by the gauge field affects world-line 
geometry. On the infinite plane, there is a simple geo- 
metrical interpretation of this interaction in terms of the 
winding numbers of the boson world lines. The wind- 
ing number w r around a plaquette r is the number of 
times the imaginary-time world lines of all the bosons 
wind around the plaquette. Consider the partition func- 
tion before averaging over the gauge field. The effect 
of the gauge field enters the partition function as the 
phase factor exp[— iJ2 a I a ' ^ x "] m ^q- (0) over tne 
gauge field. This phase factor can be written in terms 
of w r : Yl a J a • dx Q = ^2 r w r Q r . We can now perform 
the average directly over the Gaussian flux distribution 
(|Io|), instead of the gauge field distribution ([)]). We will 
be working with periodic boundary conditions (i.e., on a 
torus). This will be well-defined if we impose a constraint 
of zero total flux through the system. On averaging, the 
phase factor becomes: 



S 2 = 2tt 



-S 2 



(19) 



Thus we see that the action cost due to the current in- 
teraction is proportional to a geometrical property of the 
world lines, similar to an unoriepted area, which has been 
termed the "Amperean area"E£t 



il'r 



(20) 
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This geometrical interpretation of S2 is particularly use- 
ful in the numerical evaluation of this quantity. 

If we are working with periodic boundary conditions, 
the geometrical definition of w r given above will not work 
because there is an ambiguity in identifying which pla- 
quettes are inside or outside a loop on a torus. Neverthe- 
less, we can still use the above analysis for paths which 
do not wrap around the boundaries. (We will discuss 
wrapping paths in the next section.) The only modifica- 
tion is that we need a definition of the winding numbers 
which preserves Stokes' theorem: <f a(x) • chi = ^ r w r <fr r . 
In the case of zero total flux, a suitable definition is: 
w r = §[a®(x) — a^(x)] • dx, where a°(x) is the vec- 
tor potential at x due to a test flux $ placed at plaquette 
r, and R is an arbitrary reference plaquette. Geometri- 
cally, this picks R to be on the "outside" of any loop on 
the torus. The Amperean area as defined above is inde- 
pendent of the choice of this plaquette, because different 
choices amount to global changes in the winding numbers 
{e.g., w r — > w T + 1) and the above definition is invariant 
under such changes. 

The effect of the gauge field on the particles is now 
clear. The action S2 suppresses world-line loops with 
large winding numbers. Indeed, since S2 is non-negative, 
it excludes all configurations with finite Amperean area 
in the limit of infinite (& 2 ). This suppression can be re- 
lated to the original problem of holes moving in a spin 
liquid with a slowly varying spin quantization axis. A 
hole moving in a loop comes back with a random phase 
due to the locally fluctuating spin chiralities of the spin 
backgrounded The random phase can be interpreted as 
arising from a fictitious random flux. World-line loops 
that enclose large areas are strongly suppressed when av- 
eraged over random flux distribution due to the destruc- 
tive interference of the random phases. Therefore, we 
expect that, in the presence of strong random flux, the 
dominant contribution to the partition function comes 
from a special kind of paths that do not "see" the ran- 
dom flux, i.e., paths where fa.- dx = 0. These are "re- 
tracing paths" where each traversal of a link on the lat- 
tice is.mtraced in the opposite direction at some point in 
timecJ'ES, and such paths have zero Amperean area. 

A similar picture pf retracing paths has been studied by 
Brinkman and RiceEZi who studied a single hole in a Mott 
insulator where the spins are treated classically. Indeed, 
studies of a single particle in a strong random flux have 
yielded a density of states .nearly identical to that of the 
Brinkman-Rice problen£jEi The Brinkman-Rice model 
gives a linear- T resistivity at high temperatures (T > t) 
but a constant scattering rate of order t. Although we 
might expect this to be applicable to our model far above 
the degeneracy temperature of the bosons, this behavior 
does not extend down to the degenerate regime relevant 
to the present problem. 

At boson densities of interest here and at low temper- 
atures, Bose statistics and particle exchange are impor- 
tant; they can give rise to behavior different from the 



single-particle Brinkman-Rice result. We shall look at 
the effect of the gauge field on the quantum exchanges 
among bosons more carefully in section IV-c. For now, we 
point out that, even in the presence of strong gauge-field 
fluctuations, the bosonic nature of the particles cannot 
be ignored because the particles can form long exchange 
cycles that retraces themselves so that an individual bo- 
son does not have to retrace its own path. This is an 
important consideration at low temperatures where the 
imaginary-time paths are long allowing for a strong de- 
gree of particle exchange. Although the system can be 
highly degenerate at low temperatures, we shall now ar- 
gue that these "Brinkman-Rice bosons" remain normal 
at all finite temperatures, due to interactions with the 
fluctuating gauge fields. 



B. Destruction of superfluidity 

We will now discuss the effect of the gauge field fluctu- 
ations on the superfluidity of the Bose system. We will 
see, as in the previous section, that this can be under- 
stood in terms of the geometrical properties of the boson 
world lines. 

A neutral Bose system with short-range interaction in 
two dimensions is a superfluid below Kosterlitz-Thouless 
temperature Tkt- The onset of superfluidity at Tkt is 
caused by the binding of vortex-antivortex pairs in the 
Bose fluid so that vortex motion does not cause phase 
slips across the system. An essential ingredient of the 
existence of the superfluid phase is a long-range logarith- 
mic attraction between the vortices and the antivortices. 
A single vortex costs infinite energy in an infinite system 
E v = (irp s /m) log(L/a) where a is a short distance cutoff 
(~ vortex core radius) and p s is the superfluid density. 
Therefore single vortices cannot exist at low tempera- 
tures. Nevertheless, the proliferation of free vortices is 
possible above Tkt because this provides a gain in en- 
tropy which also scales as log L. However, in a charged 
Bose system, screening currents causes the vortex inter- 
action to be short-ranged. In our problem, the vortex 
interaction becomes exponentially weak at distances be- 
yond X P = [T/2p s t(^ 2 )} 2$ , which can be interpreted as 
a penetration depth of the Bose fluid. Now, the creat 
of a single vortex costs a finite amount of energyc. 
E v = (irps/m) log(Ap/a). This no longer compensates 
the entropic gain from vortex-antivortex unbinding, and 
so we do not expect to see a sharp phase transition of the 
Kosterlitz-Thouless type at finite temperatures. 

One might still expect that there is a crossover tem- 
perature scale below which the vortex density will be suf- 
ficiently low that the Bose system would have strong dia- 
magnetic response. A rough estimate of this temperature 
scale using a Boltzmann weight for the vortex—density 
gives a large value for this crossover temperatureEj. How- 
ever, we will see later that, in the presence of strong gauge 
fluctuations, the diamagnetic response of the bosons re- 



G 



mains small. 

To understand the suppression of superfluidity specif- 
ically in our model, we turn to the path-integral formu- 
lation of the problem with periodic boundary conditions 
(i.e., on a torus). Ceperley and Polloclcj have shown 
that superfluidity is associated with the existence of long 
world-line cycles which wrap around the torus. The su- 
perfluid density is given by 



(W 2 ) 
4/3t 



(21) 



where W x (W y ) is the number of times the boson world 
lines wrap around the torus in the x (y) direction. In 
other words, W = ^a=i Jo d,T± a /L. In the presence of 
gauge fields, superfluidity is destroyed by the same mech- 
anism that causes the Brinkman-Rice behavior: wrap- 
ping configurations pick up random phases, and should 
be suppressed by destructive interference on averaging 
over the gauge field. The number of plaquettes whose 
random fluxes contribute to the phase picked up by a 
wrapping path should increase with increasing system 
size. For a large enough system, one might expect this 
phase to be totally random. We therefore expect this sup- 
pression to be very strong. For instance, one can evaluate 
S2 for a straight-line path which wraps around the torus 
in thew-direction. To do so, we use Eq. ( [l6| ) instead of 
Eq. §L9f) because the geometrical interpretation of S2 in 
terms of winding numbers is not applicable for wrapping 
paths. We find that such a path gives 
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D(k,0) ~ 2n 2 W 2 L 2 



!A> 2 



(22) 



To compute S2 for a more general path with wrapping 
W y , one can break it down into a wrapping path with 
the same wrapping number and a non-wrapping path 
(Fig. |^) . (5*2 will consist of the contributions of the wrap- 
ping paths and non-wrapping paths separately, as well as 
a cross-term between the two paths.) We argue that S2 
diverges for all wrapping paths in the thermodynamic 
limit, and so superfluidity is destroyed at all finite tem- 
peratures. 




(b) 




FIG. 2. (a) Projection of a world line onto the xy plane 
shows a retracing path, (b) A wrapping path, (c) Decompo- 
sition of (b) into a reference path and a non-wrapping path. 

It should be noted that we are studying a gauge model 
where the uniform part of a is set to zero. We may al- 
ternatively work with a model without this gauge fixing. 



With periodic boundary conditions, this model allows an 
arbitrary Aharonov-Bohm (AB) flux through the torus. 
This flux is related to the phase of the product of RVB pa- 
rameters (^ij ) along a (Wilson) loop which wraps around 
the torus. If we average over this AB flux assuming a uni- 
form distribution, we would find that all wrapping paths 
are strictly prohibited and n s = at all temperatures 
even for samples of finite size. We will not impose such 
a drastic condition on the wrapping paths in this work. 

We can also ask whether long-range order exists in the 
Green's function for the bosons. The Green's function 
itself (&t(r)6(0)) is not gauge invariant, and would van- 
ish on averaging over different gauges. However, we can 
study the Green's function in a fixed gauge, for example 
the transverse gauge V • a = 0. In fact, one can write a 
gauge-invariant analogue correlation function which coitL- 
cides with the Green's function in the transverse gaugeEj: 



G(r) = (6 t (r)6(0)) v . a =o 

= (&t ( r )6(0)e-*/ dV /( r ')v-a(r') ) 



(23) 



where V 2 ,,f(r') = d(r' — r) — S(r'). In the path-integral 
representation, the evaluation of G involves a world line 
originating at site r and a world line terminating at site 
at the same point in imaginary time. Note that this quan- 
tity coincides with the Green's function in the Coulomb 
gauge. Consider now the phase factor ^ Q J a-dx Q picked 
up by the world lines {x a } in the evaluation of the 
Green's function G(r) in this gauge. The random flux 
$r at a distant plaquette R (with R ^> r) has a con- 
tribution of magnitude $R./i? to the vector potential at 
a point Q near and r. The sum of the contributions 
to the vector potential at Q due to the random fluxes 
at radius R from the origin is a random vector with a 
mean squared magnitude of 27ri? x (($ 2 )/_R 2 ) ~ (<J> 2 )/i?. 
This analysis is valid for all fluxes which are at a distance 
R > r. Integrating over the contribution of such fluxes, 
the variance of the magnitude of the vector potential at 
Q scales as ($ 2 ) log(L/r). Summing over all Q near 
and r, we obtain a random phase with a divergent vari- 
ance: ($ 2 )r 2 log(i/r). Thus, averaging over the distant 
fluxes for these sites, one obtains a suppression factor 
of exp[— (<fr 2 )r 2 g(r / L)] where g(x) ~ logl/x for small x. 
This can be interpreted as a binding potential for the end 
points of G(r). We therefore do not find long-range order 
in this quantity because of the destructive interference of 
the random phases due to distant fluxes. 

We will now present numerical evidence for the sup- 
pression of superfluidity below the degeneracy tempera- 
ture Tdo of the system. A measure of the degeneracy 
temperature is the Kosterlitz-Thouless temperature of 
the system at zero flux. We-make use of the observa- 
tion of Ceperley and PollockEJ that the the probability 
of bosons to participate in the multi-particle exchange 
is about \ at Kosterlitz-Thouless transition. In other 
words, the probability P\ that a boson is in an exchange 
cycle of length 1 is about \. We estimate that, for our 
lattice bosons with density rif, = 0.25 and on-site inter- 
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action U — 4t, the degeneracy temperature Too = l.li. 
(For strong on-site repulsion, Tbo is not particularly sen- 
sitive to the value of U, e.g., Tp n = OM for U = I6t.) 

We have measured, using Eq. (|2l|), the superfluid frac- 
tion n a /rib at T — i/6 with U = At for a range of flux 
variances and for systems up to 8x8 in size (Fig. |^). We 
see that the superfluid fraction decreases with increasing 
system size. In fact, the superfluid fraction as a func- 
tion of ($ 2 )L collapses onto a single curve (Fig. in- 
set), indicating that n s (L,(3, ($ 2 )) = /(T($ 2 ),/3). Since 
fix, (3) — ► as x — > oo, we see that an arbitrarily small 
random magnetic flux would destroy superfluidity in the 
thermodynamic limit. In the language of the renormal- 
ization group, this shows that the scattering by gauge 
fields is a relevant perturbation at finite temperature. 




FIG. 3. Superfluid density vs. {& 2 ) for different system 
sizes at fit — 6. Inset: a scaling plot suggests that super- 
fluidity vanishes at ( ( I >2 }c ~ 1/L. 



C. World line geometry in the normal phase 

Having established that our system remains normal at 
low temperatures, we will now examine the geometry of 
the world lines in this normal phase in the presence of 
strong gauge fluctuations. In particular, we will look at 
the effect of the gauge fields on quantum exchange and 
imaginary-time diffusion. These are mutually related: 
imaginary-time diffusion over large distances aids quan- 
tum exchange among particles and quantum exchange fa- 
cilitates imaginary-time diffusion. For example, in a dis- 
sipative model of bosons coupled to external heat bath, a 
slow logarithmic imaginary-time diffusion is expected to 
suppress quantum exchange very strongly, remitting in an 
incoherent liquid even at zero temperaturet5E£l. In our 
case, the bosons are elastically scattered by the gauge 
fields. We find that the gauge fields have less dramatic 
effects on quantum exchange and imaginary-time diffu- 
sion. 

We have shown that the world lines retrace themselves 
in the presence of random flux. One might expect that, 
compared with the case of zero flux, this would reduce 
the distance traveled by the particles in the imaginary- 
time interval (3 before their paths must return to some 



permutation of their starting positions. This should slow 
down the imaginary-time motion of the bosons as well as 
reduce the probabilities for exchange. We find that this 
is indeed the case. 

We first look at the exchange probabilities P, of a parti- 
cle participating in an exchange cycle of i bosons. As be- 
fore, we may deduce a degeneracy temperature Td from 
the probability (l-s- Pi) for a particle to be involved in 
particle exchangeO. This degeneracy temperature is re- 
duced compared to the case of zero flux. For U = it 
at quarter filling, we find that the zero-flux degener- 
acy temperature Tbo = is reduced to Td = 0.5t at 
($ 2 ) = 0.5$§. At ^-filling, it is reduced from T D0 = 0.8* 
to Td = 0.34t. A finite Td does not imply Bose condensa- 
tion at a finite temperature. Indeed, one cannot deduce 
a superfluid transition by examining the exchange proba- 
bilities. Remarkably, in the degenerate regime below Td, 
the exchange probabilities for the cases of ($ 2 ) = and 
0.5$o are nearly identical (see Table Q). In this tempera- 
ture regime, a particle is equally likely to participate in an 
exchange cycle of any size: Pi ~ P 2 — • • ■ — P/v — 1 /N, 



(a) 




FIG. 4. Schematic world-line cycles which retrace when 
projected onto the zy-plane. Solid circles denote boson po- 
sitions at t = 0. (a) Each boson retraces its own path; (b) 
Exchange cycles with more than one boson retrace their own 
paths; (c) Two exchange cycles can retrace each others paths, 
and two wrapping paths can retrace each other to give zero 
total wrapping around the boundaries. 

We can gain a qualitative understanding of the low- 
temperature exchange probabilities, by examining how 
the suppression of Amperean area by S2 affects the ge- 
ometry of the world-line configurations. When there is 
significant quantum exchange, individual bosons do not 
have to retrace their own paths in order to minimize the 
total Amperean area of the world-line configuration of all 
the bosons. Instead, one might minimize the Amperean 
area of each world-line loop formed by several bosons in 
the same exchange cycle. We find that this is not the 
entire situation at sufficiently low temperatures. Below 
Td, the different world- line loops have strong overlap. We 
find that different cycles retrace each others' paths. (See 
Fig. ||.) Thus, although the gauge fields have a drastic 
effect on the total area enclosed by all the boson world 
lines, individual world-line cycles may enclose large ar- 
eas. One might therefore expect that some aspects of 
the world-line geometry, which are insensitive to the to- 
tal area, may indeed be very similar to the case of zero 
flux. 
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The observation that individual particles do not have 
to retrace their own paths suggests that they could dif- 
fuse a greater distance than in the single-particle case. 
One should see a reduction in the kinetic energy (if) -of 
the particles compared to the Brinkman-Rice theoryEZl. 
This is indeed the case. (See Appendix for a discussion 
of the measurement of the kinetic energy.) A single parti- 
cle with retracing paths has a band edge at — 2\fit rather 
than —At. In our system, the kinetic energy per particle 
goes below the Brinkman-Rice band edge at low temper- 
atures, approaching —At roughly linearly in temperature 
(Fig. |). Thus, we see that the strong gauge fluctuations 
do not have a large effect on some aspects of the world- 
line geometry (e.g., exchange probabilities) while having 
a dramatic influence on others (e.g., superfluidity). 
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■ p/L 2 =l/6 
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FIG. 5. Kinetic energy per particle as a function of tem- 
perature. Dashed line marks the Brinkman-Rice band edge 
for the single-particle problem. (<E> 2 ) = 0.5$o an d U = At. 



Let us now examine the imaginary-time motion of 
the particles in more detail. Ideal bosons are diffusive 
in imaginary time at all temperatures, i.e., the mean- 
squared displacement of particle a is linear in imagi- 
nary time t: R 2 (t) = ([x Q (r) - x a (0)] 2 ) = Atr for 
< r < (3/2. With repulsive interactions, there is an 
increase in the effective mass of the particle, e.g., for 
U = At at quarter filling, we find t — ► t* = 0.95t. In the 
presence of random magnetic flux, the imaginary-time 
diffusion is slowed down, and the mean-squared displace- 
ment R 2 (t) is no longer linear in r at all temperatures. 
Fig. H shows our results for the (superfluid) zero-flux case 
at temperature (3t = 9 and the case of strong random 
flux at fit — 4,6,9. Since we are working with peri- 
odic boundary conditions, we have used the definition: 
R 2 (t) = ([J^\ a (r)dr] 2 } . We can see that, whereas 
R 2 (t) has significant downward curvature at j3t = 2, it 
becomes closer to diffusive behavior as the temperature 
is lowered. However, we are unable to reach the asymp- 
totic regime where the particle has traveled far on the 
scale of the interparticle spacing over a time period of 
(3/2 (see Fig. || inset). 




12 3 4 



%t 

FIG. 6. Single-particle diffusion R 2 (t) = ([x(r) -x(0)] 2 ) in 
imaginary time for < r < (3/2. Solid lines: strong random 
flux with ($ 2 ) = 0.5$o at fit = A, 6, 9. Dashed line: zero flux 
at (3t = 9. Inset: R 2 ((3/2) for zero flux (o) and ($ 2 ) = 0.5$g 
(□); dashed line marks the squared interparticle spacing. 

In order to study the long-time behavior, we can ex- 
amine the size of the world-line exchange cycles. A cycle 
where the world lines of I particles {xi, . . . ,x/} form a 
loop can be roughly regarded as a particle traveling over 
a time interval of 1(3. Thus, the possibility of exchange 
means that a world-line cycle can travel large distances 
compared with an individual boson. In a system with 
periodic boundaries, the size Ri of the cycle is defined 
by: 




For ideal bosons, Rf should equal R 2 (t = 1(3/2) at in- 
verse temperature 1(3, and therefore should scale linearly 
with I. Fig. fj] shows Rf for a 4x4 lattice with 9 particles. 
We have measured only cycles which do not have a net 
wrapping number around the periodic boundaries so that 
we do not have contributions from cycles with different 
topologies. We see that Rf is linear in I for the cases 
of zero flux and strong random flux, although the slope 
of the case with the strong random flux is reduced sub- 
stantially. This demonstrates that the imaginary-time 
motion of the bosons is diffusive at long distances. 

These results indicate that we are probing an uncon- 
ventional phase of a Bose liquid. Although the system 
remains normal, many aspects of the imaginary-time mo- 
tion of the particles in the degenerate regime resemble 
that of a neutral Bose liquid which is a superfluid in such 
temperatures. In subsequent sections, we shall study the 
physical properties of this "strange metal" and discuss 
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the relevance to the normal state of the cuprate super- 
conductors. 
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FIG. 7. Cycle sizes Rf as a function of cycle length I for a 
6x6 lattice with 9 particles. 



V. TRANSPORT AND OPTICAL 
CONDUCTIVITY 



In this section, we will present our quantum Monte 
Carlo (QMC) results on longitudinal transport for this 
strange Bose metal. To obtain the conductivity of 
the system, we measure its imaginary-time analogue 
(T a p(iui n ) in our quantum Monte Carlo simulation: 



1 



<J a p{iu n ) = I :Il a /3(iL0 n ) 



II Q/3 (*u;„) = / (,q=oW=o(0))e <u - T *-, 



(25) 
(26) 



where j q (r) = E r jrO)e lq ' r and j r (r) = J2 a 5 ( r 
x a( T ))^7^ is the gauge-invariant current (Fig. ||). 




J 0.4 



FIG. 8. Current correlation function n xx (iw n ) for a 6x6 
lattice with 9 bosons with ($ 2 ) = 0.5<3>o an d U = At. 



The imaginary-time measurements are related to the 
real-time conductivity cr(u>) = a xx {oj) by: 



1 

2~L2 



(j q=o (r)-j q=o (0)) 



'we-" T (i(w) duj 

1 - e -P" 7T ' 



(27) 



Deducing dynamical properties (such as conductivity) 
from imaginary-time data is in general an ill-posed prob- 
lem. Several approximate methods are often used in 
the context of QMC studies. A simple method, which 
has been . used in the study of the superfluid-insulator 
transitior£3Ej, is to fit a(iuj n ) to a simple functional 
form, such as the Drude form a(ioj n ) = co/(l + |w n |rt r ). 
More generally, one can use a Pade approximant to fit an 
arbitrary number of poles and zeroes: 



a(z) 



ao + aiz + . . . a,N n z 



N„ 



b + b\z 



(28) 



This approach is particularly suitable if the scattering 
rate 1 / r tr (or the position of the pole closest to the origin 
in (Eq)) is large compared to the temperature at low tem- 
peratures. This is however not the case in our problem. 
In our system, H xx (iuj n ) is nearly constant as a function 
of n for finite n even at low temperatures, suggesting that 
l/r tr is proportional to T. (Note that H xx (n = 0) = 
in the limit of strong random flux because paths which 
wrap around the torus are strongly suppressed.) 

We have calculated the conductivity by numerical ana- 
lytic coBt«H|iation using the maximum-entropy (MaxEnt) 
methodE3cll. Eq. ( p7| ) takes the form of a linear integral 
equation: 



d(r) = / K(T,w)r(u)du}, 



(29) 



where K(r, lo) is the kernel relating the imaginary-time 
data d(r) to the corresponding response function r(uj). 
In our QMC simulations, d(r) is measured at discrete 
points t 1 = I At with mean c?;. The errors for the time 
points I and m are correlated with a covariance matrix 
C/m = ((di — di)(d m — d m )). The MaxEnt method finds 
an estimate of r(u) as the function f(to) which maximizes 
the functional: 



<t>[r{u);a] = - ]pc 2 +aS, 
where x 2 is the goodness of fit 



l.m 



(A - d l )[C- 1 ] lm (D m - d m ) 



(30) 



(31) 



with Di — J dojK(ri, w)f(w), and the "entropy" S is mea- 
sured with respect to a given default model (or measure) 
m{u>): 



S = I du 



r(u>) — m[u)) — f(uj) log 



(32) 



The variable a in Eq. ( p0| ) is a regularization parameter 
controlling the competition between the smoothness and 
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the goodness of the fit, and <p is also maximized with re- 
spect to ito. We have chosen the default model m(ui) 
to be a constant in order not to build in any bias. Our 
results are not sensitive to this choice. Details of the 
MaxEnt method are given in Refs. [i(i|-[42|. 

One can check the results of the MaxEnt inversion us- 
ing relevant sum rules. In the case of conductivity, we 
have used the sum rule 



a(uj)duj 



7T (K) 
~A~1? 



(33) 



which is the lattice version of the more familiar form in 
the continuum: J Q a(uj)duj — 7ro&/2m. In our MaxEnt 
results, this sum rule is obeyed to within 3% error. In or- 
der to obtain reliable data for the MaxEnt inversion, we 
have worked with a fine discretization in imaginary time 
(iAr < 0.1). For the lowest temperatures (T < 0.4t), we 
worked at fixed Ar and f3 /L 2 . We chose f3 oc L 2 to con- 
trol the finite-size effects because of the imaginary-time 
motion of the bosons is roughly diffusive, as discussed 
above. Our results are in fact not very sensitive to this 
choice, indicating that finite-size effects are small. For 
instance, the values of resistivity at j3t = A and rib = 1/4 
for 4x4 and 6x6 are similar within statistical error. 
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FIG. 9. Optical conductivity for 6x6 lattice with 9 bosons 
at (3t =9,6,4,2,1,0.5. <$ 2 ) = 0.5$g and U = it. 

We find that a(oj) consists of a single Drude-like peak 
(Fig. ||). Since this peak exhausts the sum rule (pBSj), its 
spectral weight is proportional to —(K). This spectral 
weight has a weak temperature dependence in this tem- 
perature range because, as already discussed, the kinetic 
energy approaches — At per particle as the temperature is 
lowered. Tphis should be contrasted with the Brinkman- 
Rice resulted for non-degenerate particles (T 3> t) where 
the weight under a(cu) decreases as (-K) <~ T~ l . 

The width of <j(ui) gives a transport scattering rate 
consistent with: l/r tr = (ksT with ( — 1.8 — 2.2 
(Fig. |l0|) . This result has been obtained for two densities 
rib = 1/4 and 1/6 so that this scattering rate appears to 
be independent of density. Again this differs from the 
Brinkman-Rice result where 1 /r tr is a constant of order t 
(as one begins to see at the highest T in Fig. |l0|). The re- 
sistivity p, given by the peak height, is consistent with a 



linear temperature dependence of pe 2 /h — (l/2iriib)T/t 
for T < 2t (Fig. [l|). We estimate a statistical error of 5% 
for p by examining fluctuations due to statistical errors 
in the measurement of the current correlation function. 
There are also systematic errors due to the smoothing of 
structures. 
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FIG. 10. Scattering rate l/r tr as a function of tempera- 
ture. Solid (hollow) symbols correspond to a boson density of 
ri6=l/4(l/6). ($ 2 > = 0.53- 2 , and U = At. 
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FIG. 11. Resistivity as a function of temperature. 
Solid(hollow) symbols correspond to a boson density of 
716=1/4(1/6). ($ 2 ) = 0.5$g and U = At. 

There appears to be a systematic deviation from the 
linear-T behavior below T — 0.3t, in particular in the 
case of quarter filling. This deviation is stronger for 
p than for l/r tr . The difference can be attributed to 
the T-dependence of the Drude weight discussed above 
which should affect the resistivity but not the relaxation 
time. We speculate that the deviation from linearity at 
the lowest temperatures may indicate the approach to 
zero-temperature critical behavior. This is beyond the 
scope of this paper. 

Our resistivity agrees, to within a factor of 2, with 
Jaklic and PrelovseliEa who provided an approximate 
diagonalization of the t-J model on 4x4 lattices and 
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found a Drude peak with width 2T. They also found a 
broad background, and interpreted it with a frequency- 
dependent scattering rate t(uj). Indeed, some authors 
have interpreted the experimental optical conductiv- 
ity as possessing a power-law tail and emphasized its 
importances. This incoherent part of the conductivity 
is absent from our boson model, and may be due to in- 
elastic scattering of the bosons with the gauge field or, 
more generally, with the fermionic degrees of freedom. 



VI. MAGNETIC RESPONSE 

We now discuss the response of this degenerate Bose 
liquid to a weak external magnetic field perpendicular to 
the plane. In the absence of the random magnetic fields, 
a Bose liquid has a strong diamagnetic response as the 
temperature is lowered towards the transition to a su- 
perfluid when it develops a Meissner response. We argue 
here that the linear response of the system to a mag- 
netic field is strongly suppressed by the gauge fluctua- 
tions. Qualitatively, this can be again understood by ex- 
amining the world-line configurations. We have already 
demonstrated that the partition function is dominated 
by world-line paths which are unaffected by the internal 
gauge fields J2 a J a • dx a = for any a. These configura- 
tions are therefore also unaffected by any external mag- 
netic fields. Thus, we see that the system has a vanishing 
linear response to magnetic fields in the limit of strong 
gauge fluctuations. For the sake of completeness, we will 
now discuss more quantitatively the magnetic response 
of the system. Relevant physical quantities are the dia- 
magnetic susceptibility xb, the Hall coefficient i?H, and 
the magnetoresistance Ap/p. 

Consider first the diamagnetic susceptibility. On the 
infinite plane, in the presence of a weak external field H, 
each world-line configuration picks up an extra factor of 
cxp[-i J2 a A cxt • cfoc Q ] = exp[-i_ff„4 ] where curlA cxt = 
H, and A = J2r w r i s the oriented area of the configura- 
tion. (In this section, we will use units where $o — 2iv.) 
Expanding this in a Taylor expansion, one can write the 
partition function Z(H) as: 

Z(H) = J 2?{x} (l-iHA - l -H 2 A 2 )j e~ s ^ 



= Z(0)[l--H 2 (A 2 o ) 



(34) 



where A is the oriented area of a world-line configuration 
and (. . .) denotes an average for the system at H = 0. 
We have assumed here that the external magnetic field 
H has negligible effect on the spectrum of the gauge fluc- 
tuations. The diamagnetic susceptibility is given by: 



Xb = -5 



1 d 2 In Z 4n 2 T 



P dH 2 



(35) 



Since A a > A by definition, we can see that, when the 
gauge fluctuations are strong so that configurations with 



zero Amperean area dominate, the system has no dia- 
magnetic response, as previously suggested in Section IV 
B. 

It should be noted that, with periodic boundary condi- 
tions, the total flux penetrating the torus is quantized in 
units of the flux quantum. One should use replace (A 2 ) 
by 4(sin 2 [ J ff A/2])/ J ffo where H = $ /L 2 is the small- 
est uniform field allowed in a torus of size L. Moreover, 
as in the case for the Amperean geometrical inter- 

pretation of the phase factor J A oxt • dx is not possible for 
paths which wrap around periodic boundaries. However, 
these wrapping configurations are strongly suppressed in 
the case of strong random flux and should give negligible 
contribution to the susceptibility. 

Magnetotransport properties can be written in terms 
of the conductivity tensor cr^ g : 



Rr w a xy j(Ha 2 xx ) 
Ap/p w -Aa xx /a xx 



(36) 
(37) 



where a xx = a xx °, and Aa xx = a xx — a xx To obtain the 
conductivity tensor, we need the current-current correla- 
tor (j Q (T)/(0)) ff . Expanding again in a Taylor series in 
H, one obtains the correlator: 



jT>{x}j a jf> {l-iHA ~\H 2 A 2 ) e" s -' 



/Z>{x}(l - iHAo - ^A 2 ) e~ s «« 

(jV) - mrfAj - \a 2 ^fA 2 ) + ■ ■ 
i-±ip<^ + ... 



(38) 



Terms such as (.i x j x A n ), {j x j v A 2 ) are zero by symmetry, 
erefore, 1 

(T xy (iuJn) 



Therefore, from (|2q), we get 



H 

Aa xx (iuj n ) 
H 2 



2ni 
|^n|*o Jo 
2ir 2 



rfre— (4 =0 (r)^ =0 (0)A) 



dre 1 



[0q=oMj* = o(0K) _ (j- =0 (r)^ =0 (0))<^)]. (39) 
Since the oriented area A can be written as 



A ° = \Y. 



z ■ (jr("r) x r)dr 



1 f 

2 Jo ^' ^ q x jq ^q=° dT 



(40) 



we see that we can related this expression for the Hall 
conductivity a xy to the more familiar one involving the 
average of three currentsc3. 

Again, we see that the magnetotransport response is 
strongly suppressed by the gauge fluctuations because it 
is sensitive to the oriented area of the world-line con- 
figurations. In principle, the quantities lma xy (u:) (from 
which we can obtain ReCT^(O) from a Kramers-Kronig 
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relation) and Aa xx (u>) can be computed in a similar way 
to the calculation of the optical conductivity. However, 
these quantities are too small to measure in the regime 
of strong gauge fluctuations that we study. 

Since we have argued that the gauge field fluctuations 
are indeed strong in the cuprates at temperatures above 
the superconducting transition, it appears that our sim- 
ple boson model with a quasistatic gauge field cannot de- 
scribe quantitatively the magnetotransport in these ma- 
terials. This result is however qualitatively consistent 
with the experimental finding that these magnetotrans- 
port properties are generally suppressed from the classi- 
cal values. To obtain a quantitative prediction for these 
properties, one may attempt to restore dynamics to the 
gauge fields. If the gauge field may relax in time, then 
the boson world lines no longer have to obey the con- 
dition of strictly retracing paths. This would allow the 
world lines to enclose a finite oriented area and hence a 
finite response to external magnetic fields. However, we 
emphasize that such an approach might not represent the 
physics completely. We believe that our model illustrates 
the general point that the influence of an external field 
on the system is strongly masked by the fluctuations of 
the internal magnetic field. 



VII. DENSITY CORRELATION FUNCTION 
A. Phase separation 

Non-interacting bosons are infinitely compressible. 
They would therefore collapse into a small region of the 
system in the presence of any quenched disorder which 
has a tail of localized states in the single-particle spec- 
trum. An analogous collapse is also found in this prob- 
lem with annealed random flux-. Such an instability was 
discussed by Feigelman et allEl who have argued that 
it occurs also in the case of interacting bosons at low 
densities, leading to a hole-rich phase and a hole-absent 
phase. They further argued a long-range Coulomb repul- 
sion would be necessary to stabilize the uniform phase. 

Within the world-line picture, one can visualize the in- 
stability of the homogeneous phase in the limit of strong 
gauge fluctuations. The condition of retracing paths in 
this limit encourages the bosons to come close to each 
other so that their paths may retrace each other. This 
will allow individual boson paths to explore a larger area 
(in imaginary time), and hence lower the kinetic energy 
of the system compared to the case with each boson has 
to retrace its own path. In the absence of any repulsive 
interactions, this effect would dominate at low temper- 
atures, making the homogeneous phase unstable to col- 
lapse. 

We find that this instability towards the formation of 
dense aggregates indeed occurs in our model in the ab- 
sence of boson repulsion, although the instability is pre- 
vented by on-site repulsion, at least for the moderate 



boson densities of interest here. We have studied the in- 
stability by examining the compressibility of the system: 
k = lim g _>o K,(q) with 



/c(q) = 



1 



Nn b Jo 



dr(n q (T)n_ q (0)), 



(41) 



where H q (r) is the Fourier transform of the boson density 
at imaginary time r. Alternatively, k — limg_>o (3S(q) / n 2 
where 5(q) is the static structure factor: 



S(q) 



L 2 



(42) 




2 3 
r) (q=r|[jt/5,7t/5]) 

FIG. 12. Static structure factor (inset) and the 
q-dependent compressibility as a function of q in the (tt,tt) 
direction for different values of U. (it — 4, (<E> 2 ) 
n b = 0.25. 



0.5$g, 



In Fig. Q we show the behavior of S(q) and «(q) for 
different values of the on-site repulsion U for a 10 x 10 lat- 
tice with 25 bosons. The structure factor 5(q) as a func- 
tion of q is qualitatively different for the cases of small 
U and large U (compared to t): S(q) for q = (|, J) is 
greater than the density n& for when the on-site inter- 
action is small. We can also look at the compressibility. 
Since we work with finite systems at fixed boson num- 
ber, we will evaluate «(q) at the smallest wavevector of 
the system as an estimate of the q = behavior. We 
see that, in the presence of random magnetic flux, the 
compressibility increases with decreasing U . This can be 
interpreted as a divergence as q — + for small U , and 
hence an instability of the homogeneous phase. (This is 
also reflected in the magnitude of the fluctuations in our 
QMC results for «(q) which grows as q — > for suffi- 
ciently small U.) However, for strong on-site repulsion, 
the density correlations show no sign of an instability at 
this density. 



B. Static structure factor 

The density fluctuations in our boson model should be 
relevant to the charge fluctuations in the full t-J model. 
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It has been pointed out that the density excitations of 
the t-J model does not resemble those of a conventional 
Fermi system. We will now compare our results with 
numerical results on the full t-J model in the literature. 

The static structure factor (42) has been calculated 
by various meansErEJ. Fig. [l3] shows the static structure 
factor for_pur boson system together with that of the 
t-J modela at T = 0.25i. We see that our results are 
qualitatively similar to the t-J model, with improving 
quantitative agreement as one approaches the hard-core 
limit (see, for example, U — 16t). We should point out 
that this dependence on U should not be as strong for the 
transport properties of the system, because the particle 
currents are not directly affected by the repulsive density 
interactions. 

It is also interesting to note that the magnitude of the 
gauge field fluctuations has a relative weak effect on 5(q) 
when the on-site repulsion U is strong. However, as we 
shall see in the next section, the dynamics of the density 
excitations is strongly modified by the interaction with 
the gauge fields. 

0.25 



0.20 



3 

&5 



0.15 
0.10 
0.05 
0.00 




tj model 
Q Q U=4t, <*'>=0 

U=4t, (4> 2 )=0.5* 2 
q-q U=l6t, <* 2 >=0 
m U=16t, (3> 2 >=0.5<& 2 

12 3 4 

t| (q=r|[7c/5,7c/5]) 



FIG. 13. Static structure factor of the boson model at den- 
sity rib = 0.2 akmg the (ir,n) direction at T = 0.25t. Asterisks: 
t-J model resultQ at electron density n = 1 — nj, — 0.8 and 
t/J = 2. 



C. Dynamic structure factor 

We now look at the dynamic structure factor 5(q, ui): 
1 



S(q,«) 



L 



_> ,die^(n q (t)n_ q (0)) 



(43) 



where n q (i) is the Fourier transform of the density in 
real time. The dynamic structure factor is related to the 
imaginary-time density-density correlation function by 

— <n q (r)n_ q (0)) = J^e-^+e-^nS^^dcj. (44) 

Again, we use MaxEnt to perform the inversion of this 
integral equation. Two sum rules can be used as a check 
of the MaxEnt procedure. 



J dw(l - e- /3w )wS'(q,w) = -^(2 - cos q x - cos q y )t 

(45) 



du S{q,uj) = -n 2 b K(q) 



o 



These are lattice versions of the /-sum rule and the com- 
pressiblity sum rule. They are satisfied within 1% error 
in our results. 



<<I> /=0 

<<£ 2 >=0.5cl>- 




FIG. 14. Dynamic structure factor of the superfluid phase 
(<$ 2 ) = 0) and the normal phase ((<I> 2 ) = 0.5$o) in the (vr, it) 
direction. U = At, T — t/6 at quarter-filling. 



Fig. |IJ shows 5(q, w) for our bosons with and with- 
out the random flux. The system in the absence of ran- 
dom flux should be a superfluid at the temperature and 
densities considered here, and therefore should possess 
well-defined phonon excitations. We see sharp phonon 
peaks in the density excitation spectrum, for instance, 
at wavevector q = (f , f ) ■ These long-lived phonon 
excitations of the superfluid phase do not survive the 
coherence-breaking effect of the gauge-field interactions. 
We find only break peaks in 5(q, u>) in the presence of 
strong random flux. 

Another effect of the presence of the gauge field is a 
reduction in the bandwidth of the density excitations. 
This might be expected because the gauge-field interac- 
tion tends to increase the compressibility of the system. 
Indeed, we see that the center of the (it, it) peak is pulled 
in from 7.6t to 6.8£. 



We also see that the dynamic structure factor has 
a simple scaling with the hole density (Fig. fj"5| ): 
5(q, uj;rih) — n/ l S , °(q, oj) holds for nn = 0.1 ~ 0.3. This 
is natural in a model of degenerate bosons where the bo- 
son density is equal to the hole density. 
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FIG. 15. Scaling of S^q, lo) with boson density in the (-k, it) 
direction. The solid(dashed) lines are S(q, ui; n^/rih for 9(6) 
bosons on a 6x6 lattice, fit = 6, U = 4t, (<E> 2 ) = 0.5* 2 ,. 

We will now compare q ur . r esults with numerical re- 
sults on the full t-J modelHcil. It should be noted that, 
although we expect the electron density excitations of 
the t-J model to be dominated by its holon component, 
there is no quantitative equivalence between the struc- 
ture factors of the t-J model and our boson-only model. 
Nevertheless, we argue that the dynamic structure fac- 
tor of our model has qualitative similarities with that of 
the t-J model. For instance, the absence of sharp peaks 
in the dynamic structure factor is also found in the t-J 
model. An obvious similarity, built into our boson model 
a priori, is the lack of any structure indicating scatter- 
ing across a Fermi surface at q = 2kp. Another feature 
is the sealing of dynamic structure factor with the hole 
density^. 
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FIG. 16. Dynamic structure factor. Solid lines denote our 
Monte Carlo results for 6x6 lattice with 9 bosons at fit — 6 
with t = 0.9to. ($ 2 ) = 0.5$o and-£7 = At. Dashed lines de- 
note exact diagonalization resultaiil for 4 holes in an 18-site 
cluster with to /J = 2.5. 

We find that S(q, u) along the (ir, ir) direction agrees 
well with an exact diagonalization study of the t-J model 
at a similar hole density, as shown in Fig. [l6|. (We 
have used a moderate rescaling of the hopping energy: 
t = OMq where to is the hopping energy in the t-J 



model.) The area under q = (f , f ) peak is larger in 
our model than in the t-J model. We believe that, as in 
the case of the static structure factor, this discrepancy 
can be improved with if we use a stronger on-site repul- 
sion. However, the structure factor does not agree with 
the t-J model along the (ir, 0) direction. It might be that 
the spectrum of the holes at zero temperature is quali- 
tatively different from the simple tight-binding spectrum 
that we have assumed here. 



VIII. CONCLUSION 

In summary, we have studied a degenerate Bose sys- 
tem which remains metallic below its degeneracy tem- 
perature due to elastic scattering with random and qua- 
sistatic gauge fields. In the path-integral picture, the 
bosons retrace their paths in the limit of strong gauge 
fluctuations in order to avoid the quantum frustration 
due to the fluctuating gauge field. We have demonstrated 
that many features of these "Brinkman-Rice bosons" in- 
deed mimic the behavior of the full t-J model and the 
normal state of the cuprate superconductors. These fea- 
tures include the linear-T dependence of the longitudinal 
scattering rate and a charge excitation spectrum which 
consists of broad incoherent structures. This model itself 
has a strongly suppressed response to external magnetic 
fields, hinting that the behavior of the system as mea- 
sured in Hall and magnetoresistance experiments have 
to understood in terms of a separate mechanism. 

It would also be interesting to understand the behavior 
of the system in the zero-temperature limit. Although 
the limit of infinite gauge fluctuations (i.e., a uniform 
flux distribution on a lattice) would strictly forbid any 
world lines to wrap around periodic boundaries, one may 
consider the case of weaker gauge fluctuations in the 
zero-temperature limit and ask whether there is a critical 
value of ($ 2 ) below which the system is a superfluid at 
zero temperature. This will involve a study of the system 
at very low temperatures near a quantum critical point. 
This is beyond the scope of this paper. 
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APPENDIX A: OPERATOR AVERAGES 

In this appendix, we discuss the evaluation of operator 
averages. In the path-integral representation of the par- 
tition function, world-line configurations are sampled ac- 
cording to a distribution proportional to exp(— — S^). 
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In our Monte Carlo scheme, we discrctizc the imaginary- 
time interval (3 into M segments of length At, and we de- 
fine a world-line configuration by the boson coordinates 
at these discrete time points: {Ro, Ri, R2, ■ ■ ■ , Rm = 
P(Ro)}- (R m denotes the coordinates of the N bosons 
at m-th time slice: R n 
permutation of the coordinates of Ro-) 

The configurations are sampled according to the prob- 
ability: 

M-l 

T({R}) = - e - S2 «*» [] PAr(R m , Rm+l), (Al) 



(x^ m \ . . . , x^™- 1 ), and Rm is a 



m=0 



where Af is a normalization constant, and pat{R, R') is 
the short-time (high-temperature) density matrix in the 
absence of gauge fields. It is given by: 



PAt(R,R') = {R\e- AT ^+ H ^\R') 



(R\e 



-\ArH v -ArH° K -\ArH L 



\R>) 



(i?| e - AT <|i?')e-^ Ar( ^ (ii)+ ^ (fl ' )) (A2) 



with Hu(R) = (R\Hu\R), and 

(ij) 

Hu = ^^n^ - 1). 



(A3) 
(A4) 



The error involved in this approximation of the density 
matrix is 0(Ar 3 ). 

Measurements which depend only on particle positions 
are simple to evaluate in this path-integral representa- 
tion. The expectation value of such a measurement O is 
given by 

(O) = Tr [d(R , Ri,R 2 , ■ ■ .)V(Ro, #1, #2, ■ ■ •)] (A5) 

where 0({R}) is the measured value for the world-line 
configuration {R}, and the trace is taken over all such 
configurations. 

Averages of operators which are non-local in position 
space are more cumbersome to evaluate. An example is 
the kinetic energy: 



(ij) 



(A6) 



The gauge field aij is defined on the link between the 
neighboring sites i and j. The Peierls factor closes the 
gap in the imaginary-time loop caused by the action of 
kinetic energy operator. Inserting the operator Hk in 
the imaginary-time slice between the m = and 1, it can 
be shown that the kinetic energy can be evaluated as: 



(H K ) = Tr 



Tr 



(RolKUArlR!) 

(Ro\Uat\Ri) [W> 
(Ro\H° K e 



-AtH 



«\Ri) 



(R \e- 



-At H" 



*\Ri) 



e 2 



^AriHuiR^-HuiKRo)) 



V{{R}) 



(A7) 



where H V {KR) is defined by H V {H K \R)} = 
Hu(KR){H K \R)}, and Uat is the short-time evolution 
operator: 



Uai 



e 2 



\ArHu -ArH a K -\ArHu 



(A8) 
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TABLE I. One, two, three, and four- boson exchange prob- 
ability for various T, (& 2 ), and U at quarter-filling. 
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u 




Pi 


P 2 


Pa 


P 4 


0.5t 


At 
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0.51 


0.23 


0.13 


0.07 


Q.5t 


At 





0.20 


0.12 


0.11 


0.11 


0.25t 


At 





0.12 


0.11 


0.11 


0.11 


0.25t 


At 


0.5 


0.26 


0.16 


0.13 


0.12 


0.25t 


16t 


0.5 


0.41 


0.21 


0.13 


0.10 


O.llt 


At 





0.11 


0.11 


0.11 


0.11 


O.llt 


At 


0.5 


0.12 


0.11 


0.11 


0.11 


O.llt 


m 


0.5 


0.12 


0.11 


0.11 


0.11 
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